Proposed definitions of the correlation energy density from a Hartree-Fock starting 
point: The two-electron Moshinsky model atom as an exactly solvable model 
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In both molecular physics and condensed matter theory, deeper understanding of the correlation 
energy density e c (r) remains a high priority. By adopting Lowdin's definition of correlation energy 
as the difference between the exact and the Hartree-Fock values, here we propose two alternative 
routes to define this. One of these involves both exact and Hartree-Fock (HF) wavefunctions, while 
the second requires a coupling constant integration. As an exact analytical example of the first 
route, we treat the two-electron model atom of Moshinsky, for which both confinement potential 
and interactions are harmonic. Though the correlation energy density e c (r) is known analytically, 
we also investigate numerically its relation to the exact ground-state density in this example. 



PACS numbers: 31.15.Ew,31.25.-v,31.25.Eb 



I. INTRODUCTION 



It is true to say that one of the remaining problems 
in both molecular physics and condensed matter theory 
is to gain deeper understanding of the correlation en- 
ergy density. We have recently approached this prob- 
lem via M0ller-Plesset (MP) perturbation theory 0,0. 
By adopting Lowdin's 0] definition of the ground-state 
correlation energy as the difference between the exact 
and the Hartree-Fock (HF) values, here we shall propose 
two, formally exact, routes to the correlation energy den- 
sity. The first of these, as in Q, starts out from the 
so-called level-shift formula but in contrast to 0, 
where low-order perturbation theory is invoked, our cen- 
tral example is exactly solvable, which means that the 
MP series has been summed to all orders. This example 
is the model two-electron atom introduced by Moshinsky 
Q and it is therefore natural enough that we pose the 
two-electron atom example formally exactly in section [TXI 
immediately below. Section UTT1 presents the exact theory 
for the Moshinsky model. Comparison is made in sec- 
tion IIVI of the total kinetic energy, including correlation 
of the Moshinsky atom with that of the (non-relativistic) 
He-like series of atomic ions for large atomic number. 
An alternative route for defining the correlation energy 
density is then proposed in section [Vj which may prove 
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to come into its own in solid-state theory rather than 
molecular physics. Section I VII constitutes a summary, 
plus proposals for future studies which should be fruit- 
ful. 



II. EXACT LEVEL-SHIFT THEORY FOR THE 
GROUND STATE OF ATOMS AND MOLECULES 

Let us consider an A-body system described by the 
Hamiltonian 

H = Ha + H!, (la) 
H Q = J dT 1 dT 2 ^{v 1 )h{r 1 ,r 2 )m{r), (lb) 

ffj = ~ J dr 1 (ir2* t (ri)1' t (r2)w(ri,r 2 )*(r2)*(ri0l,c) 

where the field operators in second quantization are de- 
fined as usual by their expressions in terms of the creation 
and annihilation operators and their standard commuta- 
tion relations: 

*(r) = ^2^k(r)a k , (2a) 
k 

*t( r ) = £*£(r)4, (2b) 
fe 

[a k ,a[,] + = S k k', [at, a k '} + = 0. (2c) 

Here, v(r\,r2) is the interaction potential between the 
particles and h(ri, r 2 ) is the kernel of a one particle op- 
erator that reduces to the kinetic energy operator in the 
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usual cases. However, it may also embody the effect of 
an external potential and other effects. The creation (an- 
nihilation) operator at (af.) is assumed to create (annihi- 
late) particles in states described by the wave-functions 
\l/fc(r). Here, as is usual, r = (x, s), with x denoting the 
particle position, and s the particle spin, and k is the col- 
lective quantum number associated with the basis states 
I^Zc}. We will assume that the free Hamiltonian is diag- 
onal in the spin variable and also that the interaction is 
spin-independent, i.e. 



h(r 1 ,r 2 ) = /i(xi,x 2 )(5 Sl , 
u(rj.,r 2 ) = v(xi,x 2 ). 



(3a) 
(3b) 



Following , let us use the level-shift formula [1] , tak- 
ing as the unperturbed problem the Fock Hamiltonian 
-ffuFj with ground state energy Eq. Then, for two elec- 
trons, we can write explicitly for the correlation energy 
density e c (ri) 

, v = / ^*(ri,r 2 )[ff-g H F]$HF(ri,r 2 )rfr 2 

where ^ is the exact ground-state wave function. It 
should be noted that the definition above of the corre- 
lation energy density, Eq. (j4|), is not unique. However, 
this circumstance is not necessarily problematic, because 
any alternative definition should lead to the same to- 
tal integrated correlation energy. This is a similar situ- 
ation, and moreover also seems to be close connected, 
with the known lack of precise definitions of the energy- 
momentum tensor for general physical systems. There- 
fore, such a property should not restrict the value and 
utility of the concept, whenever it becomes possible to 
construct a theoretical scheme in which this energy den- 
sity plays a relevant role independently of its non-unique 
definition. From Eq. ^ we then have, in an obvious 
notation 

(*|$HF>e c (ri) = E J **(ri,r 2 )$HF(ri,r 2 )dr 2 

- E J **(ri,r 2 )$HF(ri,r 2 )dr 2 ,(5) 

where in reaching Eq. ([5]) from Eq. (|4|), H has been al- 
lowed to act on \I>*, and Hhf on $hf- Of course, by 
integrating both sides of Eq. ([5]) over all ri, we obtain a 
trivial identity for E — Eq, the latter quantity being sim- 
ply the total Lowdin correlation energy E c = J e c (r)dr. 

While Eq. |(SJ) is of course, valid for the (nonrelativis- 
tic) He-like series of two-electron ions with atomic num- 
ber Z, we do not presently know 'J* and E analytically. 
Therefore in section ITTTl immediately below, we turn to il- 
lustrate Eq. ((5|) analytically by appeal to the Moshinsky 
two-electron atom model |5|. 



III. EXACT USE OF THE LEVEL-SHIFT 
FORMULA FOR THE HARMONIC MOSHINSKY 
TWO-ELECTRON MODEL 

The Moshinsky model atom has confining (external) 
potential i ( |ri | 2 + |r 2 | 2 ) and particle interaction also of 
harmonic form ifc|ri —r 2 | 2 . Using coordinates R = (ri + 
r 2 )/\/2 and r = (ri — r 2 )/-\Z2, the exact ground-state 
wavefunction takes the form 



#(r, R) 



(l + 2fc)i 



■ exp 



7T2 



1. 



R 2 I exp ( --(1 + 2fe)5r 



(6) 

Less well known is the fact that the corresponding <&hf 
wavefunction entering the key expression (O for the cor- 
relation energy density has the exact form [6| 



*HF(ri,r 2 ) = 0(ri)0(r 2 ), 



where 



■ exp 



4" 



(7) 



(8) 



It is understood that Eq. (0) only provides the spatial de- 
pendence of the ground state, whose overall antisymmet- 
ric character is to be provided by the spin dependence. 
Since the spin structure of the ground state problem is 
fixed by its singlet character, we can safely employ the 
symbol r below to denote only the spatial coordinates. 

Plots of the "overlap" (^I^hf) are already available 
as functions of the particle-particle interaction strength 
k: e.g. for k = 1 the overlap is 0.94. Inserting Eqs. ([6]) 
and into the right-hand-side of Eq. (JSJ) we obtain an 
exact result for e c (r), now spherically symmetric in the 
Moshinsky atom model. The correlation energy density, 
the HF, and exact electron densities, respectively, have 
all the analytic form 



/ \ d ( a \ 3 / 2\ 6 2 - (1 - a) 2 
e c (r) = E c (-J exp(-ar 2 ), a= y — '-, 



2 (PY exp(-/3r 2 ), = 

7r / 1 + a 



(9a) 
(9b) 



/ 9HF(r) = 2(J) 2 exp(- 7 r 2 ), 7 = Vl + k, (9c) 

a = Vl + 2k, b = 1 + a + 2 7 . (9d) 

Notice that all three functions are gaussians, albeit 
with a different decay rate exponent. Figure [T] exhibits 
the fc-dependence of the latter. At k — the interac- 
tion vanishes, and all three coefficients equal one. As the 
interaction is turned on they increase, departing slowly 
from each other. Figure [H shows the reduced correlation 
energy density e c (r)/E c (full line) together with the HF 
(long dashes) and exact (short dashes) densities normal- 
ized to one, for k = 1. It is apparent from the functional 
identity and the weak parameter divergence that there 
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FIG. 1: Decay rate for the correlation energy density (a, full 
line), and the exact (/3, short dashed) and Hartree-Fock (7, 
long dashed) one-particle densities, in terms of the coupling 
parameter k (All quantities are in atomic units). 
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FIG. 2: Radial dependence of the reduced correlation en- 
ergy density (full line), exact (short dashed) and Hartree-Fock 
(long dashed) particle densities. Here, k = 1. All quantities 
are in atomic units. 



is an intimate connection between the correlation energy 
density and the exact or HF electron densities. This re- 
sult is compatible with the definition of e c (r) as propor- 
tional to the HF one particle density, as proposed in [2J. 
Exact relations for the Moshinsky model are presented 
in the Appendix. 



IV. DENSITY DEPENDENCE OF THE TOTAL 
KINETIC ENERGIES IN THE MOSHINSKY 
MODEL ATOM AND THE He-LIKE ATOMIC 
IONS AT LARGE ATOMIC NUMBER Z 



Let us illustrate in this section two examples of physi- 
cal systems in which the total kinetic energy can be ex- 
pressed as a functional of the density. We note first that 
the total kinetic energy of the Moshinsky model treated 
in section IIIII above can be expressed exactly in terms 
solely of p(r). From 0] one knows that the kinetic en- 
ergy density t now defined from the wavefunction form 



(V*) 2 is given by 
iMosh(r) = ~p(r) 
where 



3(rf-l) 2 
2 d 



2d 



1 p(r 
-log 



P(0) 



d- 1 = 2-7T 



P(0) 



. (10) 



(11) 



As noted earlier, in the Moshinsky atom the ground state 
p(r) is purely Gaussian. Notice that, in contrast to the 
He-ions at large Z, to be discussed below, there is a 
logp(r) term. Because everything is harmonic we ex- 
pect that the potential energy U will coincide with the 
kinetic term T and that both will be a half of the total 
energy E, that is 

T =[/=—. (12) 
2 ^ ; 

Thus the correlation energy density can be obtained as 



(13) 



Here, i corr (r) = t M osh(r) - e H F(r)/2, and e H F(r)/2 is de- 
fined as a density of HF energy constructed from one of 
its expressions as spatial integrals. In such a way, the 
correlation energy density is denned as e c (r) = 2i corr (r). 



1. Kinetic energy, including correlation, of the He-like 
atomic ions with large atomic number 

Following the work of Schwartz [8j on the He-like 
atomic ions with large Z (however, still non-relativistic) 
the total kinetic energy T has been obtained by Gal, 
March and Nagy [9( as 

T= A\m [^dr^-m} 2 [ P (r)d, 

(14) 

This expression for the total kinetic energy, again in- 
cluding kinetic correlation energy is evidently determined 
therefore solely by p(r) and the derivative of logp(r). 
Thus, the total kinetic energy density f(r) takes the form 

2 



t r 



P'(r) 
p(r) 



P(r) 



r=0 



p'(r) 



p{r) 



p(r). (15) 



r=0 



In closing this section we emphasize that in contrast to 
the Moshinsky atom, the total kinetic energy including 
correlation now depends on both p(r) and Vp(r). 



V. DIFFERENTIAL SECOND-ORDER 
DENSITY MATRIX FORMULA FOR 
CORRELATION ENERGY DENSITY VIA A 
COUPLING CONSTANT INTEGRATION 

The differential level shift formula (|U) for the correla- 
tion energy density e c (r) is entirely appropriate for two- 
electron systems like the Moshinsky atom. But repeated 
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volume integrations make it unwieldy for A-electron 
problems, with N > 2. Therefore, in this section we de- 
rive an alternative route via a coupling constant integra- 
tion. This leads to a formula for the Lowdin correlation 
energy density e c (r) characterized by second-order den- 
sity matrices. Since the development of efficient meth- 
ods for the evaluation of second-order density matrices is 
in rapid progess, the reduced number of integrals to be 
evaluated can represent a helpful technical advancement, 
once the expressions for the density matrices are already 
at hand. 



A. Introductory example of uniform electron liquid 



such a formula for the Lowdin correlation energy for an 
iV-electron system without the translational invariance 
of the HEL. 



B. Coupling constant formula 

Motivated by the above HEL example, we have gone 
back to the treatment by Stanton [l3| of the Lowdin cor- 
relation energy E — £7hf in terms of such a coupling con- 
stant integration. We write the Hamiltonian as 



To point the way, let us consider the homogeneous elec- 
tron liquid (HEL). As emphasized in early work, one of 
us [HI [H| has used as the 'coupling constant' the mean 
interelectronic separation r s . A variant of Hellmaim- 
Feynman theorem enables to express the ground-state 
energy for the HEL (see e.g. [l^] for a review). In par- 
ticular, the ground-state energy per electron, £7(r s ), sat- 
isfies the virial theorem [Iol | 



H(X) =H HF + XHi, 



(18) 



2T(r s ) + U(r s ) = ~r s ^ 
dr. 



(16) 



where T and U represent kinetic and potential contribu- 
tions, respectively. As shown by March and Young [11|. 
removing T = E — U from Eq. P^|) allows the resulting 
first-order differential equation for E(r s ) to be integrated 
to yield 



E(X)=-X J' 



A 2 



dX, 



(17) 



where A = r 2 plays the role of a coupling constant 
adiabatically connecting the unperturbed and the exact 
Hamiltonians, as in Hellmann-Feynman theorem. But 
it is well known that Bhf^s) = (A/rf) — (B/r s ), with 

A = | and B = ^ {^ff 2 '• Hence E(X) - £7 HF is known 
via a coupling constant integration. In HEL, U is deter- 
mined by the diagonal element of the second-order den- 
sity matrix, and hence we seek next a generalization of 



where < A < 1 is a coupling constant adiabatically con- 
necting the HF Hamiltonian to the exact Hamiltonian 
H = H(X = 1). Stanton [l3[ then generalizes Hellmann- 
Feynman's theorem to embrace the case of an iV-particle 
Hamiltonian. Beginning with the elementary identity 
E - £7 HF = £7(1) - £7(0) = J*(dE(X)/dX)dX, Eq. (12) 
in Ref. [HI reads, in obvious notation, 



E-E HF = j [<*(A)|^.|*(A)> - ($H F |^-|$HF)]dA, 

(19) 

where *(A) is the ground-state eigenvector of the exact 
Hamiltonian H(X) with eigenvalue £7(A), and <1>hf is the 
HF ground-state with energy £7hf- The scalar products 
(. . .) obviously imply an integration over the coordinates 
of N particles. The important point to stress here is that 
instead of the quantum-mechanical average (^\H\^-bf) 
entering the level shift formula, Eq. ((4|, Eq. (fT9|) involves 
'symmetric' averages like (\E , |£T|\E') and (^hfI^I^hf)- 
These, of course, are achieved at the cost of the coupling 
constant integration. But since H(X = 1) involves only 
one- and two-body operators, all the volume integrations 
but two for an TV-electron atom, molecule of cluster can 
be achieved by use of the second-order density matrix 
r(ri,r' 1 ;r 2 ,r2) defined by [U 



J 



N(N — 1) f 

r(n,ri;ra s r2) = " / **(ri,r 2 ,r 3 , . . . rj V )*(r / 1 ,r^r 3 , . . .r N )dr 3 . . . dr N . 



(20) 



Inserting Eq. ((20)) into Eq. (fT9|) we hence find 



£7- £7 



HF 



N(N 



— J dxj dr 1 dr 2 [/ l (r / 1 ,r' 2 )-/iHF(ri,4)]r(r 1 ,r' 1 ;r 2 ,r' 2 ;A)| r , =ri)r , =r2) (21) 
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where use has been made of the following results: 

dH(X) 



Hi = H — H HF (22a) 



(*hf|^|^|$hf) = ($hf|H-^hf|*hf) =0. (22b) 

In Eq. ([2Tj) . the kernel h(Y' 1: Y 2 ) associated with the 'free' Hamiltonian has been defined in Eq. (|3a|) . whereas the 
expression of the Hartree-Fock kernel /ihf (r'i , r 2 wm be given by Eq. (|2"5|) below. 

The final step, as with the level shift formula Eq. ([4]), is to drop the volume integration over r 2 and hence to 
achieve the desired result for the Lowdin correlation energy density e c (ri) as the coupling constant integration over 
second-order density matrices as 

£c(ri) = N(N- 1) I dX J ^^i^^-^HFCri,^)]^!,^;^,^;^!,,^^, (23) 



I 

The explicit form of the HF kernel /ihf(i"i, r 2 ) in Eq. (j23|) where the sum over k runs over the filled orbitals of the 
can be written as mean field problem. After explicitly writing the spin and 

spatial dependence, this expression takes the form 

/iHF(ri,r 2 ) = h(r 1 ,r 2 ) 



-5(n-r 2 ) / dr 3 v(n,r 3 ) ^ ^*(r 3 )^ fc (r 3 ) 
^ k 

-V^r^M^ni^), (24) 



h HF (x 1 ,s 1 ;x 2 ,s 1 ) = h(xi,x 2 )6 slS2 

+S SlS2 S(x 1 - x 2 ) ^2 / dx 3«(xi,x 3 ) ^*fc(x 3 ,s 3 )^ fe (x3,S3) 

s 3 =±l ^ fc 

-»(xi,x 2 ) ^ *fc( x i, si)*fc(x 2 , s 2 ), (25) 
fc 

where the single particle orbitals satisfy the HF equations 

/* 2 .HP(r 1) r 2 )* k (r 2 ) = e fc * k (r 1 ). (26) 
Now, the second-quantized HF Hamiltonian can be expressed as 

H H f = J dr 1 dr2^ Jt {r 1 )h HF (r 1 ,r 2 ) 1 ^{r2) -^2^v k 

fc 

= €ka l' ak ~ \ Vk ( 27a ) 

fc fc 

£fc = hk + v k (27b) 
h k = y rfr x rfr 2 *^(r 1 )^(r 1 ,r 2 )* fc (r 2 ), (27c) 

v k = J rfridr 2 ^(ri)^(r 2 >(ri,r 2 )(* fc (ri)^(r 2 )-^(i-i)%(r 2 )). (27d) 

k' 



It should be noted that the field operator W(r) in the single particle mean field orbitals. It should also be re- 
above expressions are now constructed in terms of the 
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called that in general the kernel /ihfM, r' 2 ) is not simply 
given by the Fock operator. It should also incorporate 
an additive constant which implements the property 

(*hf|^hf|3»hf> = (*hf|^|$hf> = E HF . (28) 

A possible optional form for i/jjp could be the one in- 
troduced in Ref. [l[ in order to propose an improvement 
of the M0ller-Plesset perturbative expansion. In this ap- 
proach the above mentioned additive constant is not re- 
quired and the second quantized version of -Hhf becomes 
a pure bilinear form in the creation and anihilation op- 
erators. 

We can, so far, only see a way to evaluate Eq. (j23|) 
wholly analytically for the two-electron Moshinsky atom. 
But we have not pressed the details of that ourselves, 
since for this simple two-electron model the level shift 
formula has overwhelming advantages over the coupling 
constant integration formula. But with recent progress 
in evaluating correlated two-body density matrices for 
systems with N > 2, we expect Eq. (|23p to rapidly be- 
come the advantageous route to employ [141 ]. Therefore, 
in App. [Bl we present an alternative method based on 
the Green function. 




FIG. 3: Reduced correlation energy density for the Moshinsky 
model as a function of the Hartree-Fock density, for a range 
of values of the latter evaluated at the origin. All quantities 
are in atomic units. 



VI. SUMMARY AND FUTURE DIRECTIONS 

The two proposals made here for the correlation energy 
density e c (ri) are embodied in Eqs. ([5]) and ([23]) . The 
first one involves knowledge of the exact many-electron 
ground-state function 'J, which of course is generally not 
available. For the Moshinsky two-electron atom model, 
however, both W and its HF counterpart are known. Our 
results show that the functional form of the correlation 
energy density, HF and exact density is the same, in sup- 
port of a simple proportionality expression. 

The second proposal has a more general character. An 
integration over a coupling constant A involving the sec- 
ond order density matrix must be accomplished. This 
alternative seems promising for solid state applications. 
However, we delay its analysis for further extensions of 
this work. The homogenous electron liquid case of the so- 
called Sawada Hamiltonian (see [15]) would appear then 
to afford a promising starting point. 

Finally, returningto the level shift, the early analytical 
work of Schwartz [8J] on He-like atomic ions with large 
Z referred to in section 4, may provide further insight 
into the use of the differential level shift formula ([4]). 
This formula, of course, is readily generalized beyond this 
two-electron example, but applications then are likely to 
involve considerable computational effort. 
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APPENDIX A: DEPENDENCE OF REDUCED 
CORRELATION ENERGY ON THE ELECTRON 
DENSITY 

This Appendix is motivated by the simple observa- 
tion that, for the Moshinsky model, where we know both 
the exact electron density p(r) and the HF counterpart 
Phf(^), we expect the difference p{r) — phf(^) to depend 
on the reduced correlation energy density e c (r)/E c . 

We have from Eqs. (|5aj) and (|5c|) in the main text, that 

p(r) = p(0)exp(-/3r 2 ), (Al) 
Phf(i-) = i o H F(0)exp(-7r 2 ). (A2) 

Hence, 

p{v) - p HF (r) = p(0)(exp(-r 2 ))^ - p H F(0)(exp(-r 2 ))^. 

(A3) 

But we also have 

e c (r) e e (0) 2 

~Er~ = -5— (exp(-r )) . (A4) 

Hence from Eq. (|A4j) we can write 

exp( - r2 »=(Si) 4 - (A5 » 



7 



Substituting this equation into Eq. (|A3[) we find 



p(r) - phf (r) = p(0) 



= (0) 



Phf(O) 



e c (r) 



^c(O), 

(A6) 

This appears to be a functional relation between 
e c (r)/e c (0) and p(r) — phf ( r ) ■ But of course, on the RHS, 
the four quantities p(0), phf(O), j/a and (3 /a all depend 
on the interaction strength k. However, one of these four 
quantities is sufficient to carry the fingerprints of k for 
the others: we single out therefore phf(O). But p(0) and 
Phf (0) are known from normalization of the densities to 
two. Hence we can write, 



P(0) 
Phf(0) 



and also 



plus 



p(r)=p(0)(exp(-r 2 )f 



Phf(i") = p H F(0)(cxp(-r 2 )) 7 
From the latter equations it follows that 



p(r)= m (0l(£ 



Phf( i ") 



.Phf(0) 

Substituting this equation in Eq. (|A6|) yields 



(A7) 



(A8) 



(A9) 



(A10) 



APPENDIX B: COUPLING CONSTANT 
INTEGRATION IN A GREEN FUNCTION 
FORMULA 



Following the discussion of Fetter and Walecka [l6| , we 
here attempt to solve the time-independent Schrodinger 
equation for an arbitrary value of A, namely 



£T(A)*(A) =£(A)tf(A), 



(Bl) 



where the wave function \P is assumed to be normalized: 
that is (^ r (A)|^ r (A)) = 1. One immediately finds then 
from Eq. JET) that 



E(X) = <*(A)|fl-(A)|*(A)>. (B2) 
Differentiating with respect to the parameter A yields 



—E(X) = (*(\)\H 1 \*(\)), 



(B3) 



where Hi = H(X = 1) — Hup, in the present study. 
Integrating this expression produces 



1 ^<*(A)[Xff 1 (A)]*(A)>. (B4) 
o A 



f Phf (r) 

W(o) 

'I 



_ Phf (r) 
Phf(0) 

ec(r) " 
ec(0) 



e c (r) 

Cc(0) 



(AH) 



Notice that the strength of the interaction does not ap- 
pear explicitly in this formula. Figure [3] shows a three di- 
mensional plot of e c (r)/e c (0), phf(f) and phf(0), where 
the latter plays the role of the interaction strength k. 



Hence again one is led to a level-shift formula, with Eq 
denoting the HF ground state energy but now via the 
coupling constant integration in Eq. (|B4|) . 

Fetter and Walecka then display the way in which 
Eq. (|B4[) . for N electrons, can be reduced to inte- 
grations over just two vectors ri and r 2 , by means 
of a time-dependent Green function, denoted now by 
G A (ri, t; r 2 , t 1 ). This result, shown in their equation 
(7.31), then reads 



E — E 



HF 



I 

±2 



— — f dvi lim lim 

A / t'—>t+ r 2 ^ri 



J dySfa-y) (ih^-hMri^)) Tr G A (y, i; r 2 , t' 



(B5) 



where in place of the kinetic energy operator, in our case the Fock operator /ihf appears. As with the level-shift 
formula in section [TTJ we now take the differential form of Eq. (|B5|) to find the second result proposed in the present 
study for the correlation energy density e c (r), namely 



e c (n) = ± l -J^ Jim J dyS(n - y) [ih^ - h^U y)) Tr G A (y, t; v 2 ,t') 



(B6) 



While Eq. JB6) is able to furnish a local definition of cor- 



relation energy density, the explicit determination of the 
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Green function G A (ri, t; r 2 , t') represents a considerable 
challenge. 

It is worth noticing that when the Green function 
is taken in the HF approximation, since the operator 
— /ihf) furnishes the inverse of this HF Green func- 
tion, it follows that the correlation energy vanishes as it 
should do, since its definition is the difference between 
the exact energy E and the HF value -Ehf- This prop- 
erty has an interesting implication. Let us assume that 
the correlation energy density defined in Eq. @ has a 
dependence on r which closely follows the one associated 
of the HF density (or the exact one ) 0. Then, from 
definition (|L36[) it follows that the HF correction to the 
propagator does not contribute at all to the correlation 
energy. Therefore, the approximate validity of the cor- 
relation energy density formula proposed in [2| means 
that the contributions to the correlation energy density 
coming from the higher order corrections to the Green 
function, should approximately follow the behavior of ei- 
ther the total or the HF density. Therefore, the validity 
of the proposals of Q for the correlation energy density, 
directly implies that the exact (or the HF) density should 
have a close relation with the exact or the HF densities. 
This conclusion arises because the higher corrections to 
the HF density turn to be approximately proportional to 
the same exact (or HF) densities as implied by the cor- 
rectness of the definitions given in Q. This property is 
supported by the analytical results of section HTT1 in which 
a close similarity between the HF electron density and the 
correlation energy density emerged. These general issues 
are expected to be analyzed in future extensions of the 
work. 

Returning briefly to the theme of the correlation en- 



ergy density e c (ri), we have written in Eq. (|B6[) a formula 
involving the Green function G x which can be brought 
into contact with a similar result for the level shift for- 
mula. Let us note that Eq. (|B6[) . as well as Eq. (UJ), can 
always be rewritten in the form: 



(n) = J C(n,r 2 )dr 2 



(B7) 



following a similar generalization for kinetic and ex- 
change energy densities, but now in HF theory, proposed 
by March and Santamaria [13, EH • Then the level shift 
(LS) formula, Eq. (@| reads 

C^r 2 )= ^*f-^f^--- r \ (B8) 
J **$HFrfri . . .dr N y ' 

The conventional exact wavefunction (WF) theory reads 



C WF (n,r 2 ) 



/ ^*H^dr 3 . . . dr N f $^ F H$ KF dr 3 ...dr 



N 



J ^*^>dr 1 ...dr N J $^ F $HFdri ...dr N 

(B9) 

While, in general, these forms of C(ri,r2) are different, 
they must all integrate to the same / e c (ri)dri. This is 
true also for the two coupling constant integration for- 
mulae derivable from Eqs. and (|B6p . 

What we stress in this Appendix is the TV-particle char- 
acter of Eqs. (|B8|) and (|B9[) . and the Green function gen- 
eralization of C(r 1; r 2 ) following from Eq. JB5|. To date, 
of course, analytical progress is restricted to the Moshin- 
sky atom. But since we have treated this fully in the 
body of the text, we shall omit further details. 
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